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In this paper, we investigate an anisotropic model of the universe by using the modified 

Xf^ ' gravity approach which involves the non-minimal y(i?)F^-type couplings of electromagnetic 

fields to gravity. After we derive field equations by a first order variational principle from the 

C^ ' lagrangian of non-minimally coupled theory, we look for a spatially flat cosmological solution 

with a large-scale magnetic field and correspondingly some kind of anisotropy. At the end 



we estimate certain values according to observations for three parameters occurring in the 



^D ■ solutions. 
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1 Introduction 

The inflation and the late-time acceleration of the universe [I] -[7] are two most exciting obser- 
vations in cosmology. Besides, in spite of that modern cosmology is based on the cosmological 
principle which assumes that the universe is homogeneous and isotropic, there are some reports 
in the literature announcing some detection of anisotropy on spatial scales where dark energy 
becomes important [8]-[10j. There exist various scenarios to explain these observations |llj-|21j. 
In this paper, we consider the modified gravity approach which involves a coupling between 
gravity and electromagnetism in Y{R)F^ form to explain them for the following reasons: 

• The non-minimal models allow us to find the solutions satisfying solar system and cosmo- 
logical tests for some coupling parameter values |22j-|24j. They are the modified dynamical 
models in which the modification of gravity is generated from the electromagnetic fields 
and the vice versa is also valid. 

• The non-minimal couplings with RF^-iy^ye which break the conformal invariance of elec- 
tromagnetic field were investigated in various aspects |25|-|36|. It is remarkable that they 
can arise from the calculation of QED one-loop vacuum polarization on a curved back- 
ground [27] and from the effective renormalization-group improved theory [37]. Then, they 
lead to electromagnetic quantum fluctuations at the inflationary stage, which lead the in- 
flation |38j-|42j. Because of the inflation at that time, the scale of the electromagnetic 
quantum fluctuations can be stretched towards outside the Hubble horizon and they give 
rise to classical fluctuations. Thus, they can be the reason of the large scale magnetic 
fields observed in clusters of galaxies [37].[33]-[i6]. 

• Such models, /(i?)-Maxwen [17]-[1S], /(G)-Maxwen t49j, /(i?)-Yang-Mills [37J and /(G)- 
Yang-Mills gravity [50] are considered to explain both inflation and the late-time acceler- 
ation of the universe. 

Although there are many other theories of modified gravity which involve a function of curvature 
scalar the non-minimal couplings with electromagnetic fields need more investigation to appear 
new insights on gravity, electromagnetism and/or electromagnetic duality. Also, it is important 
to compare consistency of the non-minimal theories with the other models according to obser- 
vations. Especially, to find cosmological solutions with magnetic fields can be important for the 
origin of cosmic magnetic fields in the universe. In this paper we consider the non-minimal cou- 
plings of gravity with electromagnetic fields by using the method of Lagrange multipliers and 
the algebra of exterior differential forms. Then we investigate some anisotropic cosmological 
solutions and discuss our results. 



2 The Theory of Y{R)F'^ Gravity 



We will derive our field equations of the non-minimal theory by a variational principle from the 
action 



I[e\u:\,F]= L= C*l, (1) 

Jm Jm 

where e" are the orthonormal basis 1-forms, coab = —^ba are the Levi-Civita connection 1-forms, 
F is the electromagnetic field 2-form and M is the four-dimensional differentiable manifold 
endowed with the metric g = rjabe"' (8) e^, rjab = diag( — \- ++). The orientation of the manifold 
is fixed by the Hodge map *1 = e" A e^ A e^ A e^. The torsion 2-forms, T", and the curvature 
2-forms, R^'bi of the spacetime are given by the Cartan-Maurer structure equations 

^" = de" + w%Ae^ (2) 

i?«^ = da;% + w"^ A uj\. (3) 
We consider the following lagrangian density 4-form [67\ [22] 

L = ^R*l-lYiR)FA*F + XaAT'' + ii/\dF, (4) 

where k^ = 8ttG is the Newton's universal gravitational constant (c = 1), i? is the curvature 
scalar, Y is the non-minimal coupling function of R, Xa and /i respectively Lagrange multiplier 
2-forms constraining torsion to zero (T" = 0) and electromagnetic fields to the homogeneous 
case {dF = 0, then we can consider F = dA). 

The non-minimal gravitational coupling of the electromagnetic field, Y{R) ^ 1, may change 
the value of the fine structure constant, i.e. the strength of the electromagnetic coupling. 
Accordingly, some observations such as radio and optical quasar absorption lines, the anisotropy 
of the cosmic microwave background radiation and big bang nucleosynthesis etc (for a recent 
review, see {5Ij ) can give some constraints on the deviation of the non-minimal electromagnetism 
from the minimal one. 

The electromagnetic field components are read from the expansion F = ^Fabe"" A e^. We use 
the shorthand notations: e" A e^ A • • • = e"^'", LaF = Fa, ibaF = Fab, i^aR°'b = Rb, i^baR"'^ = R 
where l denotes the interior product such that i^e" = 6^. The field equations are obtained via 
the independent variations of the action with respect to {e"}, {w"f,} and {F}. 

The infinitesimal variations of the total Lagrangian density L (modulo a closed form) are given 

by 

l-Se'' A R'' A *eabc + Se'' A ^: 

AKi Zi 



5L = ^Se^ AR''''A*eabc + Se'' A^Y{R){LaFA*F-FALa*F) + 6e^ ADXa 



+56" A YR{iaR'')ib{F A *F) + huab A (e^ A A*^ - e'^ A A**) 

5uab A S"'' -5F A Y{R) *F + 6XaAT'' + 6F Adfi + 6fiAdF (5) 



where Yr = -^ and the angular momentum tensor 

T."^ = ]-D[i''\YRF ^*F)]. (6) 

The Lagrange multipher 2-forms A^ are solved uniquely from the connection variation equations 

Ca A Afe - efe A Aa = 2Eaf„ (7) 

by applying the interior product operator twice 

A° = 2.bS"'' + i(^,S^^)Ae^ (8) 

We firstly substitute the Aa's into the Je" equations, then perform some simplifications and 
finally we arrive at the modified Einstein's equation 

^e!'' a *eabc + l^i^aF A*F-FAia*F) + YR{iaR'')ibiF A *F) 

+^D[i'^d{YRF^nF"''')]A*eab = 0, (9) 

while the modified Maxwell equations read 

dF = 0, d{Y * F) = 0. (10) 

3 Cosmological Solutions 

We will be considering only a magnetic field. Thus expectedly the direction, defined by it, is 
not same as the others. Accordingly, instead of an isotropic universe we naturally investigate 
cosmological model with some kind of anisotropy. The latter can be caused by electromagnetic 
field like here, or rotation, or is directly connected with nontrivial shear of a metric. Thus we 
look for exact solutions in the form of the anisotropic metric 

g = -dt^ + a{tfdx^ + b{tf{dy'^ + dz'^) (11) 

where a{t) is the polar expansion function and b{t) is the equatorial expansion function. We 
also think of the case in which magnetic fields are mainly generated rather than electric fields 
because we are interested in the generation of large-scale magnetic fields 

F = B{t)e^^ (12) 

where e^ = dt and e^ = a{t)dx, e^ = b{t)dy, e'^ = b{t)dz are the orthonormal basis 1-forms and 
B(t) represents the time-dependent x-component of magnetic field. This situation is realized if 
Y{R) increases rapidly in time during inflation |47]. This is the case here, see Fig.S.b. Besides, 
we note one point more. In the present model, the large-scale magnetic fields may be generated 



due to the breaking of the conformal invariance of the electromagnetic field through a coupling 
with the scalar curvature, Y{R)F A *F. 

As the homogenous Maxwell equation (jlOl ) gives us 

B{t) = Bo/b\t) (13) 

where Bq is a constant, the inhomogeneous Maxwell equation (jlOl i) is satisfied automatically. 
On the other hand, the modified Einstein equation ([9]) yields a set of complicated differential 
equations. In order to find some classes of exact solutions we firstly reduce them to a much 
simpler form by taking into account the following assumption 

YrB^ = -\. (14) 



K^ 



Thus we arrive at the equations 



ath-abu = 0, --l| + ^l4=0' (15) 

a b^ 2 b^ 

where subindex t denotes derivative with respect to cosmic time. The first one gives a relation 
between the polar and equatorial scale functions as 

a = aobt (16) 

where ao arbitrary constant. Thus we obtain 

2b%u - 2bHI + K^B^Ybt = . (17) 

Although generally power-law expansion or exponential expansion is considered separately, here 
we investigate a solution putting those together in a multiplicative way 

6(t) = t"exp(/3t'") (18) 

where n, m, (3 are arbitrary parameters. As comparing our results with observational data we 
use the mean expansion function 

s{t) = {ab^f^ . (19) 

Now we will try to determine n,m, 13 parameters by using the observational results of the first, 
second and third derivatives of the mean scale factor with respect to cosmic time, namely the 
present-day value of the Hubble parameter Hq, the deceleration parameter qq and the cosmic 
jerk Jq. Their definitions are given respectively by 



We can take the approximate present-day (t = to ~ 13.7 Gyr) values for the concerned quantities 
as Ho ^ 70km(s.Mpc)"^ w 7.4 x 10"^ Gyr"\ go - -0.81 and jo » 2.16 (or -0.5 < jo < 3.9) 



|52j.|53j. In accordance with observed data we estimate the observational values as n ~ 1.063, 
P w 0.0000172 and m ~ 3.053. We notice the smallness of /? compared with others. We also 
depict the six graphs of Y,R,B,H,q vs cosmic time t and Y vs curvature scalar R for these 
values of n, m, fi by setting Bq = 1 and k = 1, see Figs 1-3. 

4 Conclusion 

We have investigated consistently the early inflation, the late-time acceleration and the anisotropy 
of the universe for the non-minimally coupled gravitation theory with electromagnetic fields. Af- 
ter casting our model by a lagrangian 4- form we obtained the variational field equations. Then 
we found a class of solution under the assumption of a spatially flat anisotropic spacetime and 
a large-scale magnetic field which is generated due to the breaking of the conformal invariance 
of the electromagnetic field through its non-minimal gravitational coupling. 

Thus we point out that the non-minimal gravitational coupling of the electromagnetic field and 
the generation of magnetic fields may be a source of the inflation of the early universe and of 
the late-time cosmic speed-up. Moreover, we give certain estimations on the related parameters 
in order that our model would exhibit a behavior consistent with the current understanding 
of the observed universe. Finally we give the six related graphs from which we deduce some 
information, (i) With help of computer, R dependence of Y is obtained via Fig.l.a and Fig.l.b 
as Y{R) ~ — l/(0.2i?^'^ — i?^'^), see Fig.3.b. (ii) At the beginning the universe was too curved, 
but today its curvature is so tiny and we live in almost the de Sitter universe, Fig.l.b. (iii) As 
magnetic fleld was very dense for a short period in very far past, it is highly weak for a long time. 
Fig. 2. a. This is consistent with observational data because its observational order of magnitude 
is micro Gauss for the scales of 500 kpc [45]. (iv) The behavior of the Hubble parameter is 
similar to the magnetic fleld's behavior, Fig.2.b. (v) The transition from cosmic speed-down to 
cosmic speed-up started at time t ^ 8.2 Gyr, Fig.3.a. 
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(a) Y{t) (b) R{t) 

Figure 1: (a) The non-minimal coupling function and (b) the curvature scalar vs cosmic time t 
(Gyr) for n = 1.063, /3 = 0.0000172, m = 3.053, Bq = 1 and k = 1. 
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(a) B{t) (b) H{t) 

Figure 2: (a) The magnetic field and (b) the Hubble parameter vs cosmic time t (Gyr) for n = 1.063, 
/3 = 0.0000172, TO = 3.053, Bq = I a.nd k = 1. 
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(a) q(t) (b) Y(R) 

Figure 3: (a) The deceleration parameter and (b) tlie non-minimal function Y(K) for n = 1.063, 
/3 = 0.0000172, TO = 3.053, Bq = 1 and k = 1. 



